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Abstract 

In the framework of the httlest Higgs model with T parity, we study the Wh/Z^ pro- 
duction in association with a T-odd (anti)quark of the first two generations at the CERN 
Large Hadron Collider up to the QCD next-to-leading order. The kinematic distributions 
of final decay products and the theoretical dependence of the cross section on the factoriza- 
tion/renormalization scale are discussed. We apply three schemes in considering the QCD 
NLO contributions and find that the QCD NLO corrections by adopting the (II) and (III) 
subtraction schemes can keep the convergence of the perturbative QCD description and 
reduce the scale uncertainty of the leading order cross section. By using these two subtrac- 
tion schemes, the QCD NLO corrections to the Wh{Zh)<1- production process enhance the 
leading order cross section with a K-factor in the range of 1.00 ~ 1.43. 
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I. Introduction 



The standard model (SM) [U [2] provides a remarkably successful description of high energy 
physics at the energy scale up to 10^ GeV. Despite its tremendous success, the mechanism of 
electroweak symmetry breaking (EWSB) remains the most prominent mystery in current parti- 
cle physics, and the smallness of the Higgs boson mass cannot be protected against perturbative 
quantum corrections. The instability of the Higgs boson mass leads to the so-called "hierarchy 
problem" between the electroweak scale, mn, and the cutoff scale of the SM, A ~ 10 TeV [3]. 
However, the cutoff scale larger than 10 TeV will lead to a large radiative correction to the 
Higgs boson mass, which needs unnatural fine-tuning to give a proper EWSB scale. To solve 
the fine-tuning problem (or hierarchy problem) has become one of the main motivations to the 
construction of physics beyond the SM. Besides the supersymmetric and extra dimensions mod- 
els, the little Higgs models [3] are very attractive theories which offer an alternative approach 
to solve the "hierarchy problem" [HE], and are proposed as one kind of EWSB models without 
fine-tuning in which the Higgs boson is naturally light as a result of nonlinearly realized sym- 
metry [S]-[TU]. The simplest version of the little Higgs models is the littlest Higgs (LH) model 
[S^, which is based on an SU{5)/ S0{5) nonlinear a model [7j. However, precision electroweak 
constraints [TI] require the LH model to characterize a large value of /, so the fine-tuning be- 
tween the cutoff scale and the electroweak scale is again needed. This problem can be solved by 
the littlest Higgs model with T parity (LHT) [T^-[TC]. 

In the LHT, the SM particles are T-even and their T parity partners are T-odd. Then the SM 
gauge bosons cannot mix with the new gauge bosons, and the electroweak precision observables 
are not modified at tree-level. The loop suppression of the corrections to low energy electroweak 
observables allows the symmetry breaking scale / to be significantly lower than 1 TeV [H]. In 
order to cancel the quadratic divergence of the Higgs boson mass contributed by top loops, an 
additional T-even heavy top-quark T^ is introduced. Then a T-odd partner T_ is also required 
to implement the T parity. The CERN Large Hadron Collider (LHC) provides an opportunity 
for searching for the new particles predicted in the LHT. Many studies on the phenomenology of 
the LHT have been presented in detail [131 dZl dSl dS] • The phenomenology of these T-odd gauge 
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bosons and fermions in the LHT is very attractive. Due to the T parity conservation, the T-odd 
particles could be produced in pair at the LHC, such as (1) T-odd fermion pair (or gauge boson 
pair) production, (2) single T-odd fermion production associated with a T-odd gauge boson. 
Of all the processes in the second type at the LHC, the WhQ- associated production, where 
g_ = u^^d-,C-,S-,u^,d-,C-^S-, has the largest cross section. Therefore, it deserves special 
attention. The production rate of the pp — )■ ZhQ- + X process is of the same order as, but 
quantitatively smaller than, that of the WhQ- production at the LHC. The WhQ- and ZhQ- 
production signals could be detected by their subsequent decays g_ — >• Whq' , Wh — >• AhW and 
Zh — AhH, where q' represents u, d, c, s, n, d, c, s and Ah is the T parity partner of photon 
being undetectable as a candidate of dark matter [20J. The phenomenology of the T-odd SU{2) 
doublet particle productions at the LHC has been already studied at the leading-order (LO) 
[17| . As we know, the LO predictions for the processes at hadron colliders are always sensitive 
to the factorization and renormalization scales {fif and ^u^). For the LO Wh{Zh)q- production 
at the LHC fif enters solely through the parton distribution functions (PDFs), while the parton 
level cross section does not depend on fi^ at this order. In general, the high order contributions 
can reduce the scale uncertainty of the LO cross section, and the QCD next-to-leading order 
(NLO) corrections enhance the LO cross section. Therefore, it is important to take into account 
the QCD NLO corrections to reduce the sensitivity to these scales. 

In this paper, we focus on the Wh/Zh production in association with a T-odd (anti)quark of 
the first two generations at the LHC, pp — Wh{Zh)<1- + X (q^ = n_, c_ , s_, -U-, d_ , c_, s_), 
up to the QCD NLO. The paper is organized as follows. In Sec. II, we briefly review the relevant 
masses and couplings in the LHT. The detailed strategies of calculation are given in Sec. HI. The 
numerical results and discussions are presented in Sec. IV. Finally we present a short summary. 

II. Related LHT theory 

Based on an SU{5)/ SO{5) global symmetry breaking pattern, a subgroup [SU(2) x U{l)]i x 
[SU{2) X U{1)]2 of the SU{5) global symmetry is gauged in the LH model [15], and the gauge 
fields and Bi^ [a = 1,2,3, i = 1,2) are introduced. The kinetic terms for the gauge and 
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scalar fields can be written as 

2 



(2.1) 



where E is the nonlinear a model field of the LH model and the covariant derivative -D^S is 
defined as 



.7=1 

To implement T parity in the LHT, we make the following T parity assignment: 



(2.2) 



^IM ^ ^2m> ^IM ^ ^2/., 

n — >-^mn, where J^ = diag(l, 1,-1, 1,1). 



(2.3) 



The invariance of the above Lagrangian under T parity implies that the gauge couplings of the 
two SU{2) X U (1) gauge groups have to be equal, i.e., g\= g2 = ^/^g, 91 = 92 = V^9' ■ 

The gauge symmetry [SU (2) x U (l)]i x [SU (2) x U (1)]2 breaks down to its diagonal subgroup 
SU{2) x U{\), which is generated by the combinations {(Qf + Q|)/V2, Yi +12}, where {Qf , Yi} 
(a = 1, 2, 3, z = 1, 2) are the generators of the two SU (2) x U (1) gauge groups respectively. This 
subgroup is identified with the SM electroweak gauge group SU (2) l x C/(l)y . The corresponding 
gauge fields of the residual gauge symmetry, SU{2)l x ?7(l)y, are just the T-even eigenstates 
of the gauge sector: 



^a_Wl + Wl B, + B2 

WL-^y^ and B,-—^. 

In addition, the other four orthogonal linear combinations of gauge fields, 

Bi — B2 



H 



^ 2 and Bh 



V2 " - - ^ ' 

are odd under T parity. 

The mass eigenstates of the gauge sector in the LHT are expressed as 

Bl 



Li 



V2 



V2 



Al 
Zl 

Ah 
Zh 



cos 9w sin Ow 
— sin 9w cos 9w 

cos Oh — sin Oh 
sin Oh cos Oh 



wl 

Bh 
^3 



(T - even), 
(T - odd). 



(2.4) 



(2.5) 



(2.6) 



where the mixing angle Oh at the 0{v / f ) is defined as 

and V ~ 246 GeV is the vacuum expectation value of the SM Higgs. 

The T-even gauge bosons Al, Zl and Wl are identified with the SM photon, Z-boson and 
W^-boson, respectively. The masses of the T-odd gauge bosons are given by 

1 / 5w^\ f lv'^\ 

"lAff ^ ;^ff7 ( 1 - ) ' "iw^ff - 5/ ( 1 - g j2 ) ' mzit^rnwH- (2-8) 

To implement T parity in the fermion sector (Here we only present the description of the 
quark sector as a representative.), we introduce two incomplete SU{5) multiplets and an SO{5) 
multiplet: 

/ ^1 \ / \ / i^HR 

^-1 = , ^2=0, ^HR = XHR 

V / V V'2 / V Vhr 

A = -r\ = -r\ui, dif , {i = l,2,HR), (2.9) 

which transform under T parity as ^'i — > — So*2, *2 — > —^o^i and ^hr — > —^hr, where 
r* (i = 1,2,3) are Pauli matrices and Sq is a 5 x 5 symmetric tensor defined as 

(l2x2 \ 
1 . (2.10) 

l2x2 / 

The transformations for ^i, ^2 and ^hr under SU{5) are as ^1 — > V*^i, ^2 — > V^2 
and ^hr — > U^hr- There V G SU{5) and U is an 50(5) transformation in a nonlinear 
representation of SU (5) defined as 

C V^U^ = [/^SoF^Eo (under the SU{5) transformation V). (2.11) 

It tells us that qi, q2 and qhr are all SU{2) doublets. Therefore, we obtain two T-odd SU{2) 
doublets, qh = {qi + q2)/V2 and qHR, which are left- and right-handed respectively, and a 

T-even left-handed SU{2) doublet, qsM = (<Zi ~ Q2)/V2- 
Through the Lagrangian 

Cf = -Kf(^2^ + ^i^o^C^^)'^HR + h.c, (2.12) 
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the T-odd Dirac fermion doublet g-, defined as (g-)^ = qh and {q-)ji = Qhr, gains a mass 
equal to \/2k/ before EWSB. After EWSB, a small mass splitting between the T-odd up- and 
down-type quarks is induced, and the masses are given by [I3]-[T6] 



mu_ ^ V2Kf (^1 - , md_ = V2Kf. (2.13) 

The T-even left-handed SU{2) doublet qsM is identified with the left-handed SM fermion dou- 
blet. It can acquire Dirac masses and nid via Yukawa interactions with the T-even right- 
handed SU{2) singlets ur and dn, respectively. As we know, the up-type quark is heavier than 
the down-type quark for each generation in the SM, while the partners of the SM quarks in a 
new physics model may exhibit an inverted mass hierarchy. For example, the bottom-squarks 
are considered to be heavier than the top-squarks in the MSSM, a minimal supersymmetric 
extension of the SM. As we expected, Eq. ()2.13p indicates that the mass of the T-odd down-type 
quark is larger than that of the T-odd up- type quark in the LHT. 

In order to cancel the large quadratic divergent corrections to the Higgs boson mass induced 
by the top quark, the Yukawa interaction for the top sector must be modified. The ^'i and ^'2 
multiplets for the top sector must be completed to representations of the SU{'i)i and S'f7(3)2 
subgroups of SU{5) by introducing two additional left-handed SU{2) singlets Uli and Ul2- 
These multiplets are 

/ ^1 \ / \ 

Qi=\ Uli 1 and Q2 = ( C/l2 1 , (2.14) 

which obey the same transformation laws under T parity and SU{5) as do ^I'l and *I'2. In 
addition to the T-even right-handed SU{2) singlet uji, the top sector contains two right-handed 
SU{2) singlets Uri and Upa^ which transform under T parity as Uri < — > — C//j2- By using these 
new SU{2) singlets introduced in the top sector, we obtain four additional T parity eigenstates: 

Ul± = and Ur± = -j= . (2.15) 

Therefore, two new heavy partners with opposite T parity, T+ and T_, should appear in the top 
sector. 

As Ul- and Ur^ do not mix with uh and uhr, where uh is the up component of the 
SU{2) doublet qh, T_ is simply given by (T_)j^ = Ul- and {T^)r = Ur-. However, the T-even 



eigenstates [7^+ and Urj^ mix with usm and ur respectively, where usm is the up component 
of the SU{2) doublet qsM-, so that the mass eigenstates of the top quark t and its heavy partner 
T+ are given by 

( tL \ ^ f cos6l - sin 6'l \ / usm \ 
\ (T+)^ J V sin^L cos^L J \Ul+ y ' 

where the mixing angles 0l,_r and masses of T± are determined by the Yukawa interaction 
Lagrangian for the top sector. 

The couplings of the T-odd SU(2) doublet quarks and gauge bosons to the T-even SM 
particles used in our calculations are listed in Table [T] [131 EI], where {VHu)ij and {VHd)ij are 
the matrix elements of the CKM-like unitary mixing matrices Vhu and Vffrf, respectively. The 
two mixing matrices satisfy V^^VHd = Vckm therefore, they cannot simultaneously be set 
to the identity. In the following calculations we take Vhu to be a unit matrix, then we have 

= VcKM- 



Interaction 


Feynman rule 


Interaction 


Feynman rule 






Zi^did^ 


K-'-^-'-^){VHdhl'PL 



Table 1: The related LHT Feynman rules used in this work, where g_ = u_, (i_ , c_, s_, t_, 6_, 
i and j are the generation indices and C|f = 1 — 

III. Analytic calculations 
III..1 LO cross sections 

The partonic processes, which contribute to the parent process pp — )■ Wh{Zh)(1- + ^ at the 
LHC, are written as 

g{vi) + qi.P2) Vh{p3) + q'-iPA), {Vh = Wh, Zh, q = u,d, c, s, u, J, c, s), 

{q'_ = U-,d-,C-,S-,U-,d-,C-,S-). (3-1) 
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(1) (2) 
Figure 1: The LO Feynman diagrams for the partonic process gu — t- W^d-. 

There are two LO Feynman diagrams for each of the above partonic processes. We plot the 
LO Feynman diagrams for the partonic process gu — )• W^d- as a representative in FigHJ The 
Figmi) and Fig{T]^2) diagrams are s- and t-channel Feynman diagrams respectively. The LO 
cross section for the partonic process gq — t- VhqL has the form as 

^Lois,gq ^Ynq'-) = ——y^ y^\MLo\'^d<^2, {q = u,d,c,s,u,d,c,s), (3.2) 
4|pi|Vs J ^ 

where d^2 is the two-body phase space element, and pi is the momentum of the initial gluon 
in the center-of-mass system. The integration is performed over the two-body phase space of 
the final particles Vu and The summation is taken over the spins and colors of the initial 
and final states, and the bar over the summation indicates averaging over the intrinsic degrees 
of freedom of initial partons. 

The LO total cross section for the parent process — t- Vnq- + X can be expressed as 

CTLoiPP Vnq- + X) = 

^ <^ / dXAdXB [Gg/AixA, fJ'f)Gg/B{xB, fJ'f)^L0{9q VHq-,XAXBS,flf,^ir) + {A^ B)]\ , 



q=u,d,c,s 



(3.3) 



where Gj/p {i = g,q, P = A,B) represents the PDF of parton i in proton P, xp {P = A^B) 
is the momentum fraction of a parton (gluon or quark) in proton P, and and iir are the 
factorization and renormalization scales, respectively. 
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III.. 2 QCD NLO corrections 
III. .2.1 General description 

The QCD NLO corrections to the pp — )• Wh{Zh)q- + X process involve the following compo- 
nents: 

(i) The QCD one-loop virtual corrections to the partonic processes gq — )• Wh{Zh)(i'-- 

(ii) The contributions of the real gluon emission partonic processes gq — )• WH{ZH)qL + g- 

(iii) The contributions of the real light-quark emission partonic processes gg — ?■ WH{ZH)q'-+q, 
q"q" ^ WH{ZH)q'- + q and qq" ^ WH{ZH)q'- + q" ■ 

(iv) The corresponding contributions of the PDF counterterms. 

It should be noticed that for the gg — )• WH{ZH)q'- + q and q"q" WH{ZH)q'^ + q light- 
quark emission partonic processes there exists resonance effect due to the g_ propagator. In 
Figj2] we present the Feynman diagrams for these real light-quark emission partonic processes 
via intermediate on-shell T-odd quarks. To deal with the resonance effect in these partonic 
processes, we replace mg_ in the denominator of the q- propagator by m^^ — imq_Tq_. With 
the LHT parameter values used in this paper, the main decay channels of q- are q- — )• Wnq' , 
q- — > Znq and q- — )• Anq- 

Br(g_ ^ Wnq) + Br(g_ ^ Znq) + Br(g_ ^ Anq) ^ 100%. (3.4) 

Therefore, the value of Tq_ is obtained approximately by summing up the LO partial decay 
widths of these main decay channels. These QCD NLO contribution parts from the gg — )• 
WH{ZH)q'- + q and q"q" — )• WH{ZH)qL + q partonic processes are quite large due to the high 
gluon luminosity and the q- resonance effect. 

In this work, we apply three schemes in considering the QCD NLO corrections. In scheme 
(I) (denoted as "QCD NLO I") we include all the four components mentioned above in the 
QCD NLO corrections. With this scheme, the Feynman diagrams in FigjJ] could lead to large 
corrections to the Born pp — )• WH{ZH)q-+X process due to the q- resonance effect, and destroy 
the perturbative convergence. Furthermore, these Feynman diagrams are also counted towards 

9 
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Q-Q- production followed by an on-shell decay — t- Wh{Zh)q' ■ Therefore, to avoid double 
counting and to keep the convergence of the perturbative QCD description of the Wh{Zh)q- 
associated production channel, we should remove the intermediate on-shell g_ contributions 
from the Wh{Zh)q- associated production [23]. 

In scheme (II) (denoted as "QCD NLO 11") we exclude the contributions of the partonic 
processes gg — WH{ZH)qL+q and q"q" — )• Wh{Zh)<i'-+(1 from the QCD NLO corrections. Since 
the corrections to the parent process pp — t- Wh{Zh)q- +X contributed by the partonic processes 
gg —7- Wh{Zh)q'- + Q and q"q" — t- Wh{Zh)q'- + Q and their corresponding PDF counterterms 
are IR-safe, we could exclude them from the QCD NLO corrections to the parent process pp — t- 
Wh{Zh)(1- + X. With this subtraction scheme, the intermediate on-shell g_ contributions to 
the Wh{Zh)q~ associated production are removed and the perturbative convergence is kept. 
Since all the gg and qq initiated contributions are excluded, this scheme subtracts some genuine 
QCD NLO contributions. 

We adopt another subtraction strategy, the PROSPINO scheme \22\ 123). which removes the 
on-shell T-odd quark pair production from the real light-quark emissions gg — )• Wh{Zh)q'^ + Q 
and q"q" — )■ Wh{Zh)q'- + Q, to avoid double counting and to not artificially ruin the convergence 
of the perturbative QCD description of the pp — >• Wh{Zh)<1- +X process. This on-shell subtrac- 
tion scheme can provide a reliable production rate since it only subtracts the squared on-shell 
amplitudes and does this point by point over the entire phase space. We call this subtraction 
scheme as the scheme (III) in this work, which is defined as a replacement of the Breit-Wigner 
propagator [23] 

\M\\sv,,) ^ .35) 

(sy^g -m2J2 + m2_r2_ (^svhq - + ml_Vl_ 

\M\'^{ml_) 2 

where sy^g is the squared momentum flowing through the intermediate q_ propagator. The 
results by adopting this scheme are denoted as " QCD NLO III" . 

To isolate the UV and IR singularities, we adopt the dimensional regularization method 
in D = 4 — 2e dimensions. The collinear counterterm of the PDF, 6Gi/p{x, fif) {P = proton, 
i = g,u, u, d, d, c, c, s, s), is split into two parts: the collinear gluon emission part 6G^f^p°"'\x, ^if) 
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9 dn du 
(3) 



Figure 3: The box Fynman diagrams for the partonic process gu — t- W^d- 



and the coUinear hght-quark emission part dG 



^Gq{g)/p{x,fif) = (^gJ^^"°p^(x,^/) + {q = u,u,d,d,c,c,s,s), (3.6) 



(quark) , 
i/P ^ 



(quark) , 



where 



6G^ 



5G^^-'\x,^,f) = \ 



as r(l - e) ( 47r/i^ 
2^r(l-2e) [~Jif^ 

as r(l -e) / ATTfif 
27rr(l -2e) /i^ 

r(l - e) / 47r;u2 



1 



Pgg(99)(z)Gg(^gyp{x/z,Hf), 



1 



Pqg{z)Gglp{x/z,Hf), 



d,d,c,CyS,s I 



/ -Pgq{z)Gq/p{x/z,^lf). (3.7) 



q=u,u 



27rr(l - 2e) I fij 

The exphcit expressions for the sphtting functions Pij{z) {ij = qq,qg, gq, gg) in Eqs. (l3.7p are 
available in Ref. 



III. .2. 2 Virtual and real emission corrections to gq — Wh{Zh)<i'- 

The one-loop level amplitudes for the partonic processes gq — )• WHiZH)q'- in the LHT include 
the contributions of the self-energy, vertex and box graphs. In Figl3]the box Feynman diagrams 
for the partonic process gu — )• W^d- are presented as a representative. 

The strong coupling constant, the masses and wave functions of the relevant colored particles 
in the LHT are renormalized to remove the UV divergences of the virtual corrections. In the 
QCD NLO calculations of the pp — t- Wh{Zh)<1- + ^ process, the following renormalization 
constants are introduced: 



0,L,R 
«(<?-) 



l+\sz^ G„ 



ni 



mq_ + 6mq_ , 



9s = 9s + S9s, 



(3.8) 
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where Qs denotes the strong couphng constant, niq^ is the T-odd quark mass, Vg^'^ ) and 
denote the fields of the SM quark, T-odd heavy quark and gluon, respectively. The masses and 
wave functions of the colored fields are renormalized by adopting the on-shell scheme, then the 
relevant renormalization constants are expressed as 



6Z, 



5Z„ 



5Z, 



Sit 

Qs(/»r) 



Auv - ^IR 
Ac/y + 2A/R + 4 + 31n 



j4_ 



Srun 



5Za 



Sir 
27r 



Auv + In 



+ 4 



(3.9) 
(3.10) 

(3.11) 



(g_ = u^,d^,C-, S-,t-,b-), (3.12) 



where Auv = — 1e + ln(47r) and Ajji = l/ejR — je + ln(47r). 



For the renormalization of the strong coupling constant gs , we adopt the MS scheme at the 
renormalization scale ^r, except that the divergences associated with the massive top-quark, T- 
odd SU (2) doublet quarks {u- ,d-,C-,S-,t^,b^) and T± loops are subtracted at zero momentum 
]. Then the renormalization constant of the strong coupling constant can be obtained as 



9s 



«s(^r 



T- 



3 ^ 1 , mr 1 , mi, 1 x - , 

-Ac/y + -ln^ + - 2^ ln^ + -2^1n 



mi 



T=r+ ' ' q- 
{q- = U-,d-, C-, S-,t-,b-). 



The LO amplitude for gq — )• Wh{Zh)(i'- can be expressed as 



Mlo = Ms + Mu 



(3.13) 



(3.14) 



where Ms and Mt are the amplitudes for the s- and i-channel Feynman diagrams, respectively. 
Then the QCD NLO counterterm amplitude can be written as 



^as , 1 

gs 2' 



McT = ( — + -5Zy + -5Zg + -5Zq,_ Mlo + bm^_Mt 



1,-, 
2 



(3.15) 



where M.i 



represents the amplitude obtained from the t-channel amplitude 



-m / y 



Ait by doing the replacement of 



f'q'_ -'^q'_ ) 



, and p„' is the four-momentum of the 
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T-odd quark in the f-channel propagator. From Eqs. (j3.12p and (j3.13p we can see that the terms 

T- ^2 

of Yl 1^ ~f exactly canceled in Eq. p.l5p . therefore, the values of rriT^ are unnecessary in 

T=T+ 

our numerical calculations. 

We use our developed in-house programs to isolate analytically the IR singularities of loop 
integrals and calculate numerically one-loop integrals based on the LoopTools-2.4 package \26\ 
[27] , where the analytical expressions for the IR-singular parts of loop integrals are adopted from 
Ref.[28j, and the numerical evaluations of IR-safe A^-point {N < 4) integrals are implemented 
by using the formulas in Refs. [29l [30t |3T]. 

We employ the two cutoff phase space slicing (TCPSS) method [23] to calculate the correc- 
tions from the real gluon/light-quark emission partonic processes. An arbitrary soft cutoff 6s 
separates the real gluon emission phase space into two regions, the soft gluon region and the 
hard gluon region. Another cutoff 6c decomposes the real hard gluon/light-quark emission phase 
space region into the hard collinear (HC) region and the hard noncollinear {HC) region. Then 
the soft and collinear IR singularities are isolated from the IR-safe region. The integration over 
the HC region of phase space is performed in the four-dimensions by using the Monte Carlo 
integrator |32j . Finally, the total cross section for the real emission process can be expressed as 

Ao-i? = Ao-s Ao-H = Ao-s + Ao-HC + ActjJc^. (3.16) 

The UV singularities of the loop corrections are canceled by those of the related counterterms 
contributed by the renormalization constants in Eqs. (|3.9p - (|3.13p . Therefore, the renormalized 
virtual corrections (loop corrections combined with the related counterterms) are UV-finite. 
Furthermore, the renormalized virtual corrections also contain soft and collinear IR singularities. 
These IR singularities exactly vanish after combining the renormalized virtual corrections with 
the contributions of the real gluon/light-quark emission processes and the PDF counterterms. 
These cancelations have been verified analytically and numerically in our calculations. 



14 



IV. Numerical results and discussions 
IV.. 1 Input parameters 

In the study of the dependence of the QCD NLO corrected cross section on the factorization 
and renormahzation scales, we set the two unphysical scales equal to a common value (a*/ = 
Hr = /u) and do not vary them in an independent way for simplicity. This setting of scales 
may render the results more stable than they actually are due to the logarithmic term In ^ 
in the QCD NLO contributions. For example, the PDF counterterms in Eqs. (j3.7p have the 
form as ^( - — 7^; + ln47r + In ^ ) (P tS) G), where P G represents the convolution of the 
splitting function P with the PDF G. When we set fif = /i^, we obtain In ^ = and the 
factorization/renormalization scale dependence of these PDF counterterms is underestimated. 

We take one-loop and two-loop running in the LO and QCD NLO calculations, respectively 
[33j . The central value of the factorization/renormalization scale /i is chosen as = {rriWif + 
We adopt the CTEQ6L1 and CTEQ6M parton densities with five flavors in the LO 
and NLO calculations, respectively [34] . The strong coupling constant as{^) is determined by 
the QCD parameter A^'^ = 165 MeV for the CTEQ6L1 at the LO and = 226 MeV 

for the CTEQ6M at the NLO [33]. We ignore the masses of u-, d-, c-, s-, 6-quarks, and take 
aewiml)-^\j^ = 127.925, mw = 80.399 GeV, mz = 91.1876 GeV, mt = 171.2 GeV and 
sin2 0^ = l- (^)' = 0.222646. 

The colliding energy in the proton-proton center-of-mass system is taken as y/s = 7 TeV for 
the early LHC and y/s = 14 TeV for the later running at the LHC. The Cabibbo-Kobayashi- 

Maskawa (CKM) matrix elements are taken as 

/ Kd Vus Vub \ I 0.97418 0.22577 \ 
VcKM = V,d Vcs Vcb = -0.22577 0.97418 . (4.1) 

V Vtd Vs Vtb J \ 1 / 

IV.. 2 Independence on two cutoffs 

In order to demonstrate the independence of the total QCD NLO correction to the pp ^ ug ^ 
W^d- + X process on the two cutoffs, 6s and 6c, we present the QCD NLO correction parts 
as the functions of the cutoffs in FiglH^a), where we take 6c = 6s /WO, f = 600 GeV and 
K = 1. From Eqs. (|2.8| ) and (j2.13p . we obtain the related masses of the T-odd particles as 

15 



= 398.57 GeV, mu_ = m^^ = 830.70 GeV and m^, = ms_ = 848.53 GeV . In this figure 
we take fi = fj,Q = {rriwH + rnd_)/2 = 623.55 GeV. Although the decay width of is less than 
1% of mq_, we take Tii_ = O.lm^. to suppress the resonance effect, which makes the cutoff 
independence more clear. The amplified curve for Aatot of FigUJ^a) is shown in Figl^b). The 
figures demonstrate that the total QCD NLO correction Aatot which is the summation of the 
two-body and three-body corrections, is independent of the two cutoffs within the statistical 
errors, even though the two-body correction (Acj^^^) and three-body correction {Aa^^'^) are 
strongly influenced by the cutoffs 6s and 6c- As we know, the independence of the total QCD 
NLO correction to the pp ^ ug ^ W^d-. + X process on the cutoffs 6s and 6c is a necessary 
condition that must be fulfilled for the correctness of our calculations. In the further numerical 
calculations, we Rx 6s = 1 x 10^^ and 6c = 6s /WO. 



b 

<l 









tot 




Vs=14TeV,S =5/100 

c s 




f=600GeV,K=1 Aa^* ..... 









-6.0 -5.5 -5.0 -4.5 

Log(5) 



-4.0 



0.130 

b 
<1 

0,170 



(b) 



t i t i I 1 



Vs=14TeV,S =8/100 

' c s 

f=600GeV,K=1 



-6.0 -5.5 -5.0 -4.5 

Log(5) 



-4.0 



Figure 4: (a) The dependence of the QCD NLO corrections to the pp ^ ug ^ W^^d^ + X 
process on the cutoffs 5s and 6c at the LHC, where / = 600 GeV , k = 1, 6c = (^s/100 and 
fj, = fj,Q = {m\Yjj + "irf_)/2 = 623.55 GeV . (b) The amplified curve for Aatot of FigUJ^a). 



IV.. 3 Dependence on factorization/renormalization scale 

In Figsi5]^a,b,c) and Figsl6]^a,b,c) we present the LO, QCD NLO corrected cross sections and 
the corresponding K-factors for the pp — )• WhQ- + X and pp — )■ Z^q^ + X processes as the 
functions of the factorization/renormalization scale at the LHC with ^/s = 7 TeV and 14 TeV , 
respectively. In FigsEUa,b) and FigsE]^a,b) the LHT input parameters are taken as / = 500 GeV 
and K = 1, while in FigEl^c) and FigEJc) we take / = 1 TeV and k = 1. The masses of Wh, Zh 
and q- {q = u, d, c, s) corresponding to these LHT parameters are presented in Tabled In these 
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K 


/ 




rriu^ = mc_ 


'm.d_ = ms_ 






(GeV) 


{GeV) 


{GeV) 


{GeV) 


{GeV) 




500 


322.1 


685.7 


707.1 


514.6 




700 


457.8 


974.7 


989.9 


723.9 


1 


900 


592.3 


1260.9 


1272.8 


932.5 




1000 


659.3 


1403.5 


1414.2 


1036.7 




1100 


726.1 


1545.9 


1555.6 


1140.9 




1300 


859.7 


1830.3 


1838.5 


1349.1 




500 


322.1 


2057.1 


2121.3 


1221.7 


3 


700 


457.8 


2924.0 


2969.9 


1713.8 




900 


592.3 


3782.7 


3818.4 


2205.3 



Table 2: The masses of Wh, Zh and g_ {q = u, d, c, s) for some typical values of the LHT 
parameters. 

figures the curves labeled by "NLO I", "NLO 11" and "NLO III" are for the QCD NLO corrected 
cross sections using the (I), (II) and (III) schemes, respectively. The figures show that by using 
the (II) and (III) subtraction schemes we can get almost the same and moderate QCD NLO 
corrections to the production rate with a strongly reduced factorization/renormalization scale 
uncertainty in the plotted range of /i, while the QCD NLO corrections using the scheme (I) do 
not obviously improve the scale dependence of the LO cross section and destroy the perturbative 
convergence in some range of //. In the following analysis we set the factorization/renormalization 
scale /i as its central value /io = {^Wh + "^(i_)/2- 

IV.. 4 Dependence on LHT parameters 

We depict the LO, QCD NLO corrected cross sections and the corresponding K-factors for the 
pp — )• Wnq^+X a.ndpp — )• Zhq^+X processes as the functions of /, the SU{5) global symmetry 
breaking scale of the LHT, at the LHC with y/s = 7 TeV and 14 TeV in Figsi7|^a,b,c) and 
Figsj8{a,b,c), respectively. In Figsi71^a,b) and Figsl8]^a,b) the parameter k is set to be 1, while in 
FigiD^c) and FiglH^c) we take k = 3. The curves labeled by "NLO I", "NLO 11" and "NLO III" 
are for the QCD NLO corrected cross sections using the (I), (II) and (III) schemes, respectively. 
One can conclude from these figures that the cross section for the pp — )• Wh{Zh)q~ + X process 
decreases quickly with the increment of /, because the two final T-odd particles become heavier 
with the increment of /. However, in the plotted range of / we could have observable production 
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Figure 5: The dependence of the cross sections and the corresponding K-factors for the pp — >■ 
WhQ- + X process on the factorization/renormahzation scale fj, at the LHC. (a) / = 500 GeV, 
K = 1 and i/i = 7 TeV. (b) / = 500 GeV, k = 1 and = 14 TeV. (c) / = 1 TeV, k = 1 and 
= 14 TeV. 
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Figure 6: The dependence of the cross sections and the corresponding K-factors for the pp — >■ 
ZhQ- +X process on the factorization/renormahzation scale /i at the LHC. (a) / = 500 GeV, 
K = 1 and i/s = 7 TeV. (b) / = 500 GeV, k = 1 and ^/s = U TeV. (c) / = 1 TeV, k = 1 and 
= 14 TeV. 
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rates for the pp WhQ- + X and pp — )• Z}iq_ + X processes, especially when k = 1. 



566 707 849 990 1131 1273 1414 566 707 849 990 1131 1273 1414 1556 1697 1838 1980 2121 




566 707 849 990 1131 1273 



m,(GeV) 




f(GeV) 1 .0 i • 1 • 1 • 1 • 1 

' ' 400 500 600 700 800 900 

m„ (GeV) . . ^ ^ ^ . 

" 253 322 390 458 525 592 



Figure 7: The cross sections and the corresponding K-factors for the — )• WhQ- + X process 
as the functions of the LHT parameter / at the LHC. The corresponding rriyy^ and m^^ values 
are also scaled on the x-axis. (a) k = 1 and ^/s = 7 TeV. (b) k = 1 and ^/s = 14 TeV. (c) 
K = 3 and y/s = 14 TeV. 



We present the numerical results for the pp — )• WhQ- + X and pp — t- ZhQ-^ + X processes at 
the LHC for some typical values of the LHT parameters in Tableland Tabled respectively. In 
the two tables we list the values of the proton-proton colliding energy y^, the LHT parameters 
K and /, the LO cross section ctlo^ the QCD NLO corrected cross sections using the (I), (H) and 
(ni) schemes, and the corresponding K-factors. From both the tables and the /-dependence 
figures, we can see that the QCD NLO corrected cross sections for the pp — )• WhQ- + X and 
pp — )• ZhQ-^ + X processes at the LHC by using the subtraction scheme (II) are almost the same 
as those by adopting the subtraction scheme (HI) in the LHT parameter space considered in 
this paper. At the early LHC, the QCD NLO corrected cross section for the pp — )• WhQ- + X 
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Figure 8: The cross sections and the corresponding K-factors for the pp — )• ZhQ- + X process 
as the functions of the LHT parameter / at the LHC. The corresponding m^y^ and md_ values 
are also scaled on the x-axis. (a) «; = 1 and y/s = 7 TeV . (b) k = 1 and y/s = 14 TeV . (c) 
K = 3 and ^ = 14 TeV. 
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K 


/ 


0"LO 


"nlo 


Kin 


"nlo 




"nlo 




(TeV) 




[GeV) 




ifb) 




ifb) 




ifb) 








500 


253.53(1) 


544.6(8) 


2.15 


296.6(8) 


1.17 


295.3(8) 


1.16 


7 


1 


700 


21.721(1) 


40.52(7) 


1.87 


27.87(7) 


1.28 


27.73(7) 


1.28 






900 


2.5287(1) 


4.352(9) 


1.72 


3.514(9) 


1.39 


3.503(9) 


1.39 






500 


2830.7(1) 


7648(3) 


2.70 


2823(2) 


1.00 


2832(3) 


1.00 






700 


432.10(2) 


1011.9(4) 


2.34 


467.9(3) 


1.08 


465.8(3) 


1.08 


14 


1 


900 


93.171(3) 


197.66(8) 


2.12 


106.86(7) 


1.14 


106.45(8) 


1.14 






1100 


24.460(1) 


47.82(2) 


1.96 


29.32(2) 


1.20 


29.08(2) 


1.19 






1300 


7.2751(3) 


13.349(6) 


1.83 


9.053(6) 


1.24 


8.995(7) 


1.24 






500 


26.543(1) 


46.14(6) 


1.74 


31.42(5) 


1.18 


31.42(6) 


1.18 


14 


3 


700 


1.5903(1) 


2.320(6) 


1.46 


2.072(6) 


1.30 


2.067(6) 


1.30 






900 


0.12741(1) 


0.1841(5) 


1.44 


0.1810(5) 


1.42 


0.1806(5) 


1.42 



Table 3: The numerical results for the pp — >• WhQ- + X process at the LHC for some 
typical values of the LHT parameters. 



{pp —7- ZhQ- +X) process using the subtraction scheme (III) can reach 295.3/6 (158.2/6) and the 
corresponding K-factor is 1.16 (1.20) when / = 500 GeV and k = 1. While at the 14 TeV LHC, 
the QCD NLO corrected cross section for the pp — WhQ- + X {pp — ZhQ- + X) process using 
the subtraction scheme (III) can reach 2832 fb (1491/6) and the corresponding K-factor is 1.00 
(1.03) when / = 500 GeV and k = 1. We can also find that the QCD NLO K-factor increases 
with the increment of the LHT parameters n and /. By adopting the subtraction scheme (HI), 
the K-factor can reach 1.42 and 1.43 when k = 3 and / = 900 GeV for the pp — )• WhQ- + X 
and pp — )■ ZhQ- + X processes at the 14 TeV LHC, respectively. 

IV.. 5 Transverse momentum and rapidity distributions of final particles 

In this subsection we inspect the characteristics of the transverse momentum and rapidity distri- 
butions of the final decay products. The WhQ- associated production at the LHC can be followed 
by the subsequent decays Wh AhW, q- Wrq' AnWq' and Therefore, 
the WhQ- production signal can be found by detecting the final states + jet + -ET.missing 
(-ET^missing = A.hAhvi')- Similarly, the ZhQ- production can be detected through the de- 
cays Zh — AhH, H ^ bh and g_ — )• Whq' AnWq' — )• Anli^q', with the final states as 

/^66 + jet + i?T,missing (-£^T,missing = Ah Ah V ). 

The results in subsections IIV..3I and IIV..4I show that the QCD NLO corrections using the 
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K 


/ 




"nlo 


kW 


"nlo 




"nlo 




(TeV) 




[GeV) 


ifh) 


ifb) 




ifb) 




ifb) 








500 


132.15(1) 


282.7(6) 


2.14 


158.9(6) 


1.20 


158.2(7) 


1.20 


7 


1 


700 


11.205(1) 


20.94(5) 


1.87 


14.54(5) 


1.30 


14.54(6) 


1.30 






900 


1.2977(1) 


2.232(6) 


1.72 


1.813(6) 


1.40 


1.808(6) 


1.39 






500 


1446.5(1) 


3898(6) 


2.69 


1497(6) 


1.03 


1491(6) 


1.03 






700 


219.36(2) 


513.6(9) 


2.34 


241.6(9) 


1.10 


240.6(9) 


1.10 


14 


1 


900 


47.141(4) 


99.9(2) 


2.12 


54.6(2) 


1.16 


54.3(2) 


1.15 






1100 


12.351(1) 


24.22(6) 


1.96 


14.89(6) 


1.21 


14.83(6) 


1.20 






1300 


3.6685(3) 


6.758(9) 


1.84 


4.598(9) 


1.25 


4.582(9) 


1.25 






500 


14.125(1) 


24.54(7) 


1.74 


17.02(6) 


1.21 


17.18(7) 


1.22 


14 


3 


700 


0.83522(7) 


1.227(4) 


1.47 


1.095(5) 


1.31 


1.102(4) 


1.32 






900 


0.066474(5) 


0.0969(3) 


1.46 


0.0947(3) 


1.42 


0.0952(3) 


1.43 



Table 4: The numerical results for the pp — >• ZhQ- + X process at the LHC for some 
typical values of the LHT parameters. 



scheme (I) would destroy the convergence of the perturbative QCD description of the Wh{Zh)<1- 



associated production, while the QCD corrections using the (II) and (III) subtraction schemes 



are almost the same and can keep the perturbative convergence. In the following discussions on 



the transverse momentum {px) and rapidity (y) distributions of final particles, we present only 



the numerical results by adopting the subtraction scheme (III). 



The LO and QCD NLO corrected transverse momentum distributions of the final W - 



boson, jet, lepton / and missing energy {AhAhvi') for the pp — )• WhQ- + X process and the 



corresponding K-factors are plotted in Figsl9]^a,b,c,d), respectively. The results in these figures 



are obtained by taking / = 1 TeV, n = 1 and ^/s = 14 TeV. With these parameters we obtain 



niAif ~ 153.0 GeV from Eq. (j2.8p . It should be declared that in Figl9t|b) the px distribution 



labeled by "NLO III" is for the leading jet {ji is called as leading jet if Ej-^ > Ej^.), if there 



exist two iets in one event. From these figures we find that the distributions j — - — , -j—^ — and 
— — increase with the increment of px in the low px region, and reach their maxima at 

'-''PT, missing 

Px,w- ~ 200 GeV, pTjet ~ 500 GeV and 

PT,missing ~ 400 GeV , respectively. The transverse 



momentum distribution of / is quite different from the former ones. It decreases rapidly with 



the increment of Pxi- for experimentally acceptable lepton. 



The LO and QCD NLO corrected transverse momentum distributions of the final H^-boson, 
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0.08 




200 400 600 800 200 400 600 800 1000 



Figure 9: The LO and QCD NLO corrected pT distributions of final particles for the pp — > 
Wuq- + X process at the LHC by taking / = 1 TeF, k = 1 and a/s = 14 TeV. (a) W-boson, 
(b) jet, (c) lepton (d) missing py. 
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jet, lepton and missing energy {Ah Ah v ) for thepp — )■ process and the corresponding 

K- factors are depicted in FigsITOlfa.b.c.d). respectively. We do not distinguish the electric charge 
of and in the "PT distributions of the final VF-boson and lepton. In these figures we set 
the input parameters as / = 1 TeV ^ k = 1 and y/s = 14 TeV. Again the px distribution labeled 
by "NLO III" in FigllOlfb) is for the leading jet in a two-jet event. These figures show that the 
transverse momentum distributions of W^, jet, and missing energy for the pp — )• ZhQ- + X 
process are quite similar to those of W~ , jet, and missing energy for the pp — )■ WhQ- + X 
process, respectively, while the production rate of the pp — )• WhQ- + X process is almost twice 
larger than that of the pp — )• ZhQ- + X process. 



0.04 




200 400 600 800 200 400 600 300 1000 



PT,,(2eV) p^^,,,,„^(GeV) 

Figure 10: The LO and QCD NLO corrected px distributions of final particles for the pp — )• 
ZhQ- + X process at the LHC by taking / = 1 TeV, k = 1 and y/s = 14 TeV. (a) VF-boson, 
(b) jet, (c) lepton, (d) missing pT- 



We present the LO and QCD NLO corrected rapidity distributions of final particles and the 
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corresponding K-factors for the pp — )• Wh<1-+X and pp — )• Wnq^+X processes in FigsJTTTa.b.c) 



and FigsJTWa.b.c). respectively. Because 



da 
dy 



da 
dy 



for the processes at the LHC, we plot 



only the distributions . The values of the input parameters are taken the same as those used 
in FigsiH] and FigsJlOl For the ZhQ- production, we do not distinguish the electric charge of 
and in the rapidity distributions of the final VF-boson and lepton. We can see from these 
figures that the rapidity distributions ^^^"^ | , ^[^^ and ^|^°" | for the pp — )• WhQ- + X process 
are quite similar to the distributions ^j^^j d\y^ t\ ^^"^ ^^Tf ^'^^ ~^ ^hQ- + X process, 

respectively. All these differential cross sections decrease with the increment of \y\. That means 
all the final products, including VF-boson, jet and lepton, prefer producing transversely. 





Figure 11: The LO and QCD NLO corrected rapidity distributions of final particles for the 
pp — )• WhQ- + X process at the LHC by taking / = 1 TeV, k = 1 and ^/s = 14 TeV. (a) 
VF~-boson, (b) jet, (c) lepton . 
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Figure 12: The LO and QCD NLO corrected rapidity distributions of final particles for the 
pp — > ZhQ- + X process at the LHC by taking / = 1 TeF, k = 1 and ^/s = 14 TeV. (a) 
IF-boson, (b) jet, (c) lepton. 
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V. Summary 

In this paper, we calculate the Wh{Zh)(1- {q- = u-,d-,c-, s-,u-,d-,c-,s-) associated pro- 
duction at the CERN LHC up to the QCD NLO. We investigate theoretically the dependence 
of the cross section on the factorization/renormalization scale, and present the transverse mo- 
mentum and rapidity distributions of the final decay products. 

The real light-quark emission partonic processes gg WH{ZH)q-+q and q"q" — )• WH{ZH)q-+ 
q are at the same order of the on-shell T-odd quark pair production with the subsequent decay 
gr_ Wh{Zh) + q- Including these partonic processes will give a large contribution to the NLO 
QCD corrected rate for the associated WH{ZH)q- production. In order to show how to provide a 
reliable production rate of the pp WH{ZH)q- +X process at the LHC, we adopt three schemes 
in considering the QCD NLO corrections for comparison. Our calculations demonstrate that by 
using the scheme (I) the perturbative convergence could be destroyed due to the double counting 
showing up, while we can keep the convergence of the perturbative QCD description and obtain 
moderate QCD NLO corrections to the production rate with a strongly reduced scale uncertainty 
by adopting the (II) and (III) subtraction schemes. The smallness of the discrepancy between 
the numerical results using the (II) and (III) subtraction schemes indicates that the contribu- 
tions of the partonic processes gg — > WH{ZH)q- + q and q"q" — > WH{ZH)q'- + q are dominated 
by their g_ resonance effects. The QCD NLO corrections by adopting these two subtraction 
schemes enhance the LO cross section with a K-factor in the range of 1.00 ~ 1.43. We conclude 
that for associated production processes like the pp — > Wnq- + X and pp — > Znq- + X processes 
investigated in this paper, it is crucial to implement a consistent and reliable on-shell subtraction 
scheme separating associated production from QCD mediated pair production properly. The 
scheme (II) subtracts some genuine QCD NLO contributions, since all the gg and qq initiated 
contributions are excluded from the QCD NLO corrections. Therefore, The PROSPINO scheme 
is more consistent that the subtraction scheme (II). 
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